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Abstract. 

The transition from hadronic matter to quark matter in the core of neutron stars is likely to be associated with the appearance 
of a mixed phase, leading to a smooth variation of the star density profile. We discuss the results of a systematic study of the 
properties of the mixed pase upon Coulomb and surface effects. A state of the art nonrelativistic equation of state of nuclear 
matter has been used for the low density phase, while quark matter has been described within the MIT bag model, including 
the effect of perturbative one-gluon exchange interactions. The implications for neutron star structure are discussed. 
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1. Introduction 

The possible occurrence of a core of deconfined quark matter 
in the interior of neutron stars has been investigated by a num- 
ber of authors over the past three decades (for a recent review 
see, e.g., Weber |2004i. Due to the complexity of the under- 
lying dynamics, theoretical approaches to the study of hybrid 
stars largely rely on models to describe the equation of state 
(EOS) of strongly interacting matter, in both the hadronic and 
quark sector, as well as on a set of assumptions on the nature 
of the phase transition. 

Most calculations have been carried out using nuclear EOS 
obtained from either nonrelativistic nuclear many body the- 
' ory (NMBT) or relativistic mean field theory (RMFT), while 
the deconfined phase is usually described within the MIT bag 
model (Chodos et al [l974l . 

In their pioneering work Baym and Chin ( 1962 1 employed 
the familiar Maxwell double tangent construction (see, e.g., 
Huang 1963 1, which amounts to assume that the transition oc- 
cur at constant pressure. Within this picture, charge neutral nu- 
clear matter at energy density e^M cohexists with charge neu- 
tral quark matter at energy density £qm, the two phases being 
separated by a sharp interface. 

In the early 90s Glendenning (Glendenning [19921 [19971 
first pointed out that the requirement that the two phases be 
individually charge neutral is in fact too restrictive. In a more 
general scenario charged nuclear and quark matter may share a 
common lepton background, thus giving rise to a mixed phase 
extending in space over a sizable fraction of the star. 

The appearance of a mixed phase strongly affects the 
macroscopic properties of the star. A transition at constant pres- 
sure necessarily leads to the appearance of a discontinuity in 
the density profile, i.e. to a star consisting of a inner core of 



quark matter at energy density eqm surrounded by nuclear mat- 
ter at energy density 6nm- On the other hand, the mixed phase 
allows for a smooth variation of the energy density, leading in 
turn to a smooth variation of the star density profile. 

Whether the transition proceeds at constant pressure or ac- 
cording to Glendenning 's picture depends upon i) the value of 
the Debye screening length, driving charge separation, and ii) 
the amount of electrostatic and surface energy needed for the 
formation of the structures of quark and nuclear matter in the 
mixed phase. If screening is too strong to allow for a uniform 
lepton density, or if the energy loss due to Coulomb and surface 
effects exceeds the gain in bulk energy, the standard scenario 
predicted by Maxwell construction turns out to be favorable. 

The results of a detailed study carried out by Heiselberg et 
al ( 1993 1, suggest that the mixed phase is energetically favored 
over a sizable density range if the surface tension is less that 
~ 70 MeV/fm 2 . In their calculations these authors adopted the 
MIT bag model (Chodos et al ll974i for the quark matter EOS, 
while nuclear matter was described using a somewhat oversim- 
plified model, in which the energy-density included a quadratic 
compressional term and a symmetry term taken from a previ- 
ous work of Lattimer et al ( 199"Tl. 

Over the past decade, the availability of new nucleon- 
nucleon potentials, resulting from accurate fits to nucleon- 
nucleon scattering data, and the improvement of the compu- 
tational schemes made it possible to obtain a new generation of 
EOS within the framework of NMBT ( Akmal & Pandharipande 
[19971 Akmal et al [T998l . In view of the fact that NMBT is a pa- 
rameter free approach, whose dynamics is strongly cosntrained 
by nuclear data, and has been shown to possess a highly re- 
markable predictive power in theoretical studies of few nucleon 
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systems (Pieper & Wiringa 1995 1, it provides a natural candi- 
date to describe neutron star matter in the nuclear phase. 

In this paper we extend the work of Heiselberg et al ( 1993 1 
carrying out a systematic study of the stabilty of the mixed 
phase. We adopt the state of the art EOS of nuclear matter, 
obtained by Akmal et al ( 1998 1 within NMBT, and analyze the 
dependence of the results on i) the parameters entering the MIT 
bag model EOS, employed to describe quark matter, and ii) the 
value of the surface tension, driving both Coulomb and surface 
effects. 

In Section|5]we summarize the main features of the model 
EOS of both nuclear and quark matter, while the implementa- 
tion of Gibbs conditions in the case of two chemical potentials, 
leading to the appearance of the mixed phase, is discussed in 
Section|3] Our main results are presented in Section|4l devoted 
to the role of Coulomb and surface energy. The implication of 
the appearance of the mixed phase for neutron star structure are 
outlined in Section [5] Finally, the conclusion of our work are 
stated in Section|6] 

2. Models of neutron star matter EOS 

In this Section we summarize the main features of the EOS 
employed in our work, focusing on the region of nuclear and 
supranuclear density (riB ^ 0.1 frrT 3 ). For the lower density 
region, corresponding to the outer and inner crust of the star, 
we have used the EOS of Baym et al ( 1971 ) and Pethick et 
al {T995 1, respectively. It has to be pointed out, however, that 
our results are largely unaffected by the details of the EOS at 
subnuclear density, as the fraction of star mass in the crust is 
only about 2 %. 

The theoretical descriptions of nuclear and quark matter are 
both based on the standard assumptions that the system be at 
zero temperature and transparent to neutrinos produced in weak 
interaction processes (see, e.g., Shapiro & Teukolski ll997> . 

2.1. Nuclear matter 

Within NMBT, nuclear matter is viewed as a collection of 
pointlike protons and neutrons, whose dynamics is described 
by the hamiltonian 



Zj ?M Zj 1 



j» 



k>j>i 



(1) 



M and p,- being the nucleon mass and the momentum of the i-th 
particle, respectively. The nucleon-nucleon (NN) potential v; 7 
reduces to the well known Yukawa one-pion exchange potential 
at large distances, while its short and intermediate range behav- 
ior is determined by fitting the available experimental data on 
the two-nucleon system (deuteron properties and ~ 4000 NN 
scattering phase shitfs) (Wiringa et al 1995 1. The three-nucleon 
potential Vy*, whose inclusion is needed to reproduce the ob- 
served binding energies of the three-nucleon system, consists 
of the Fujita-Miyazawa two-pion exchange potential supple- 
mented by a purely phenomenological repulsive interaction ( 
Pudliner et al H995l . 



The many-body Schrodinger equation associated with the 
hamiltonian of Eq.([0 can be solved exactly, using stochastic 
methods, for nuclei of mass number A < 10. The resulting ener- 
gies of the ground and low-lying excited states are in excellent 
agreement with experimental data (Wiringa & Pieper 1995 1. 
Exploiting translational invariance, accurate calculations can 
also be carried out for uniform nuclear matter (Wiringa et al 
IT9881 Akmal & Pandharipande [T997l . 

Akmal & Pandharipande ( 1997 > have used cluster expan- 
sions and chain summation techiniques to obtain the energy 
per particle of both pure neutron matter (PNM) and symmetric 
nuclear matter (SNM). In their approach the Argonne vi$ po- 
tential of Wiringa et al ( 1995 1 is modified to take into account 
the fact that NN potentials fitted to scattering data describe in- 
teractions between nucleons in their center of mass frame, in 
which the total momentum P vanishes. Within the approach of 
Akmal & Pandharipande relativistic corrections arising from 
the boost to a frame in which P + 0, are included up to order 
P 2 /M 2 . 

Interpolating between the PNM and SNM results of Akmal 
& Pandharipande ( 1997 1, Akmal et al ( 1998 1 have determined 
the energy of matter with arbitrary ptoton fraction, needed to 
obtain the EOS of yS-stable matter, consisting of neutrons, pro- 
tons electrons and muons. Their calculations span a range of 
baryon number density iib extending up to ~ 8 no, «o = 0.16 
fm -3 being the empirical saturation density of symmetric nu- 
clear matter. 

At any given value of ns proton and lepton densities are 
determined by the requirements of charge neutrality, 



and /3-stability, 

He = H„ ■ 



(2) 



(3) 



(4) 



In the above equations n,- and yu, denote the number density 
and chemical potential of the particle of type i (i - n,p,e,fi), 
respectively. 

The EOS P - P(e), where P and e denote pressure and 
energy density is obtained from the density dependence of the 
binding energy per baryon, E = E(ns), through 



P = T 



n 



1 dE 

2 dn 



B 



and 



e = n B (E + M) 



(5) 



(6) 



The results discussed in the following sections have been 
obtained using the EOS of Akmal et al (1998), hereafter re- 
ferred to as APR, to describe the nuclear phase of neutron star 
matter. The corresponding energy per baryon of PNM, SNM 
and /3-stable matter is displayed in Fig.^ as a function of den- 
sity. 
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Fig. 1. Energy per baryon of nuclear matter calculated by 
Akmal et al ( 1998 i,plotted as a function of baryon number den- 
sity. The dashed, dotdash and solid lines correspond to pure 
neutron matter, symmetrical nuclear metter and /3-stable mat- 
ter, respectively. The box represents the equilibrium properties 
on symmetrical nuclear matter obtained from extrapolation of 
empirical data. 

2.2. Quark matter 

Due to the complexity of the fundamental theory of strong in- 
teractions (Quantum Chromo-Dynamics, or QCD) a first prin- 
ciple description of the EOS of quark matter at large density 
and low temperature is out of reach of the existing compu- 
tational approaches. To describe the quark matter phase, we 
have used the simple MIT bag model (Chodos et al l 19741 1. in 
which the main features of QCD are implemented through the 
assumptions that: i) quarks occur in color neutral clusters con- 
fined to a finite region of space (the bag), whose volume is 
limited by the pressure of the QCD vacuum (the bag constant 
B), and ii) residual interactions between quarks are weak, and 
can be treated in low order perturbation theory. 

Within the MIT bag model the thermodynamic potential £2 
can be written as 



Q = Q„ erl + VB, 



(7) 



where £2 pe n and VB, V being the normalization volume, denote 
the perturbative and nonperturbative contribution, respectively, 
and 



pert 



v 2 1 a? 



(8) 



In the above equation, the index / specifies the quark flavor, 
while Q?* is the n-th order term of the perturbative expansion 
in powers of the strong coupling constant, a s . 

The EOS of quark matter can be obtained from the relations 
linking pressure and energy density to Q: 

P = _5 = _ jB _y (9) 
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rif and pj being the density and chemical potential of the 
quarks of flavor /. 



The lowest order perturbative contributions to the thermo- 
dynamic potential, corresponding to Q^' with n — and 1, read 
(see, e.g., Tamagaki & Tatsumi 1993i 
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The chemical potentials appearing in Eas.l ll 1> and d!2t can be 
written 



fx f = e Pf + Sfif = yjmj + p 2 Ff + Sn f , 



(13) 



where the first term is the Fermi energy of a gas of noninter- 
acting quarks of mass mj at density ttf = pp /n 2 , whereas the 
second term is a perturbative correction of order a s , whose ex- 
plicit expression is ( Baym & Chin ! 19761 



6n f = 



2a s 
3^ 



Pf, - 3 — - In 



e Ff + p Ff 
m f 



(14) 



In the case of massless quarks Eqs.(|9} and (II Oi reduce to the 
simple form 



3 

4V 2 



(15) 



(16) 



(17) 



reminiscent of the EOS of noninteracting quarks. 

For any baryon desity, quark densities are dictated by the 
requirements of baryon number conservation, charge neutrality 
and weak equilibrium. In the case of two flavors, in which only 
the light up and down quarks are present, we have 



implying the relationship 



n B = \(n u + n d ) , 



— n„ 

3 



3 Hd ' 



p d =p u + p e 



(18) 

(19) 
(20) 



where n e and \i e denote the density and chemical potential of 
the electrons produced through d — * u + e~ + v e . Note that we 
have not taken into account the possible appearance of muons, 
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as in the relevant density region fi e never exceeds the muon 
mass. 

As the baryon density increases, the af-quark chemical po- 
tential reaches the value fij = m s , m s being the mass of the 
strange quark. The energy of quark matter can then be lowered 
turning d-quarks into s-quarks through d + u — > u + s. In pres- 
ence of three flavors, Eas.l ll8M20> become 



800 



"b = ^(n u + n d + n s ) , 

2 1 1 

- -n<i - -« s - n e = 

fl d = fl s = fl u + H e . 



(21) 

(22) 
(23) 



Unfortunately, the parameters entering the bag model EOS 
are only loosely constrained by phenomenology and their 
choice is somewhat arbitrary. 

As quarks are confined and not observables as individual 
particles, their masses are not directly measurable and must 
be inferred from hadron properties. The Particle Data Group 
(Hagiwara et al 2000 ) reports masses of few MeV for up and 
down quarks and 60 to 170 MeV for the strange quark. We 
have set m u = mj = and m, = 150 MeV for the up, down 
and strange quark, respectively. In the density region relevant 
to our work heavier quarks do not play any role. 

The strong coupling constant a s can be obtained from the 
renormalization group equation, yielding 

«. = 2 - . (24) 

(33 -2^/)ln(7i 2 /A 2 ) 

where Nf — 3 is the number of active flavors, A is the QCD 
scale parameter and Jl is an energy scale typical of the rele- 
vant density region (e.g. the average quark chemical potential). 
Using A ~ 100 v 200 MeV (see, e.g., Ellis et al [T996T l and set- 
ting Jl = fid ~ fi u at a typical baryon density ng ~ 4«o one gets 
a s ~ 0.4 -r 0.6. The results discussed in this paper have been 
obtained with a s = 0.5. 

The values of the bag constant resulting from fits of the 
hadron spectrum range between ~ 57 MeV/fm 3 , with A = 220 
MeV, (De Grand et al [T975l a nd ~ 350 MeV/fm 3 , with A = 
172 MeV (Carlson, et al !1983> . However, the requirement that 
the deconfinement transition do not occur at density ~ «o con- 
strains B to be larger than ~ 120- 150 MeV/fm 3 , and lattice 
results suggest a value of ~ 210 MeV/fm 3 (Satz 1982). In or- 
der to gauge the dependence of the results upon the value of B, 
we have carried out our calculations setting B = 120 and 200 
MeV/fm 3 . 

Fig-Elshows the energy density of neutral quark matter in 
weak equilibrium as a function of baryon density for differ- 
ent values of B and a s . The solid and dashed lines have been 
obtained setting a s = 0.5 and B = 200 and 120 MeV/fm 3 , 
respectively, while the dashdot line corresponds to a s = 
and B = 200 MeV/fm 3 . Comparison between the dotdash line 
and those corresponding to a s + suggests that, contrary to 
what stated by many authors (see, e.g., Steiner et al 2000), per- 
turbative gluon exchange, whose inclusion produces a sizable 
change of slope, cannot be simulated by adjusting the value of 
the bag constant and must be explicitely taken into account. 
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Fig. 2. Energy density of neutral quark matter in weak equi- 
librium as a function of baryon number density. The solid and 
dashed lines have been obtained setting a s = 0.5 and B = 200 
and 120 MeV/fm 3 , respectively, while the dashdot line corre- 
sponds to a s = and B = 200 MeV/fm 3 . The quark masses are 
m u = mj = 0, m s = 150 MeV. 



3. From nuclear matter to quark matter 

Early works on the possible occurrence of quark matter in 
neutron stars (e.g. Baym and Chin 1962) were based on the 
assumption that nuclear and quark matter were both charge 
neutral. As a consequence, the transition was described using 
Maxwell construction (see, e.g., Huang 1963 1 and the resulting 
picture of the star consisted of a quark matter core surrounded 
by a mantle of nuclear matter, the two phases being separated 
by a sharp interface. 

In the 90s, Glendenning (119921 IT9971 pointed out that this 
assumption is too restrictive. More generally, the transition can 
proceed through the formation of a mixed phase of charged 
nuclear and quark matter, global neutrality being guaranteed 
by a uniform lepton background. 

Equilibrium between charged phases of nuclear matter 
(NM) and quark matter (QM) at T = requires the fulfillment 
of Gibbs conditions (see, e.g., HuangQ%3} 



B _ B Q _ Q 

(*NM ~ ^QM ' L*NM ~ ^QM 



(25) 



(26) 



where P denotes the pressure, while fi B and ffi are the chemical 
potentials associated with the two conserved quantities, namely 
baryonic number B and electric charge Q. 

The above equations imply that, for any pressure P, the 
projection of the surfaces /Vm(jU B ,/U G ) and Pqm(i-i b ,I-i q ) onto 
the P — P plane defines two curves, whose intersection cor- 
responds to the equilibrium values of the chemical potentials. 
As the chemical potentials determine the charge densities of 
the two phases, the volume fraction occupied by quark matter, 
X, can then be obtained exploiting the requirement of global 
neutrality 



xQqm + (1-x)Q. 



NM 



2> 



. 



(27) 
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where Qqm, Qnm and Q[ denote the electric charge carried by 
nuclear matter, quark matter and leptons, respectively. From 
Eq.Cl it follows that 



X 



Qnm + Hf Qc 
Qnm - Qqm 



(28) 



with < ^ < 1. Finally, the total energy density e can be 
calculated using 



e =X^QM + (1 ~x)£nm , 



(29) 



and the EOS of state of the mixed phase can be cast in the 
standard form P - P(e). 

Requiring that the two phases be individually neutral, as 
in the pioneering work of Baym & Chin ( 1962 1, reduces the 
number of chemical potentials to one, thus leading to the equi- 
librium conditions 



Pnm(Mnm^ ~ Pqm(.^qm) 



f*NM - ^QM ■ 



(30) 



(31) 



Within this scenario, charge neutral nuclear matter at baryon 
number density n^ M coexists with charge neutral quark matter 
at density n® 
ments 



L ' M rig M and n® M being determined by the require- 



_ ( de NM \ _ l de QM \ 

^ \ dn B / ns= „r I dtl B In 



(32) 



At n N B M 



< ris < rig 



pressure and chemical potential remain 
constant, the energy density is given by 



(33) 



and the volume fraction occupied by quark matter grows lin- 
early with density according to 



X = 



^NM(n NM ) 



eQM(n® M ) - e NM {n" B M ) 



NM\ 



(34) 



Note that the above equation obviously implies that < ^ < 1, 
with^(4 M ) = 0and^(«f M )= 1. 

In the present work the intersection between the surfaces 
describing the pressure of nuclear and quark matter has been 
determined numerically, choosing as independent variables, in- 
stead of n B and yft, the proton and neutron chemical poten- 
tials /j. p and In nuclear matter they are simply related to 
the lepton chemical potential through the /3-stability condition 
(*n ~ Hp — Me- In quark matter the chemical potentials of up 
and down quarks can be obtained from fi p and /j„ inverting the 
relations 



Hp = 2n u + jj d 



(35) 
(36) 



while the strange quark and lepton chemical potentials are dic- 
tated by the conditions of weak equilibrium 



(37) 
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Fig. 3. Isobars P(ji n ,fi p ) = 200 MeV/fm 3 (upper panel) and 
350 MeV/fm 3 (lower panel) obtained using the APR EOS of 
nuclear matter (solid lines) and the MIT bag model of quark 
matter, with a s = 0.5 and B = 200 MeV/fm 3 (dashed lines). 
The intersections determine the values of the chemical poten- 
tials corresponding to equilibrium of the two phases according 
to Gibbs rules. 



(38) 



Fig- ID illustrates the construction we have employed to de- 
termine the values of fi p and /u„ corresponding to equilibrium 
between nuclear matter, described by the APR EOS, and quark 
matter, described by by the MIT bag model EOS with a s = 0.5 
and B = 200 MeV/fm 3 . The region P mi „ < P < P max in which 
the isobars of nuclear and quark matter intersect defines the 
range of densities n„„„ < n# < n max in which the mixed phase 
is energetically favored. At «g < n m ,„ and tib > n max the ground 
state consists of pure nuclear and quark matter, respectively. 

The phase transition between nuclear and quark matter, ob- 
tained setting B = 200 MeV/fm 3 and a s = 0.5, is illustrated 
in Fig.|4] Dashed and dotdash lines show the dependence upon 
riB of the energy density of charge neutral nuclear and quark 
matter in weak equilibrium, respectively, while the solid line 
corresponds to the mixed phase. The latter turns out to be the 
ground state of neutron star matter at densities .7 < iib ^ 1.7 
fm" 3 . 

The dependence of our results upon the MIT bag model 
parameters can be gauged from the upper panel of Fig. [5] It 
clearly appears that a lower value of the bag constant, corre- 
sponding to a softer quark matter EOS, leads to the appearance 
of the mixed phase at lower density. Keeping a s = 0.5 and 
setting B = 120 MeV/fm 3 one finds that the mixed phase is en- 
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Fig. 4. Dashed and dotdash lines show the energy density of 
charge neutral nuclear and quark matter in weak equilibrium, 
respectively. The bag model parameters have been set to B = 
200 MeV/fm 3 and a s = 0.5. The solid line corresponds to the 
mixed phase, obtained from Gibbs equilibrium conditions. 



ergetically favored in the range .6 < n# < 1.4 firT3. An even 
larger effect, illustrated by the lower panel of Fig.[5]is obtained 
with B = 120 MeV/fm 3 and a s = 0, i.e. neglecting perturba- 
tive gluon exchange altogether. For this case we also show the 
results obtained from Maxwell construction, leading to the co- 
hexistence of charge neutral nuclear matter at rig = .42 ftrT 3 
and charge neutral quark matter at rig = .57 fm~3. This cohex- 
istence region is to be compared to the region of stability of the 
mixed phase, corresponding to .22 < < .75 firT 3 . 



4. Structure and stability of the mixed phase 

The results of the previous Section suggest that, irrespective of 
the details of the EOS, the transition from nuclear to quark mat- 
ter proceed through the formation of a mixed phase. However, 
two issues relevant to both the appearance and the stability of 
the mixed phase, not taken into account in the discussion of 
Section|3 need to be further analyzed. 

Consider a mixed phase consisting of droplets of quark 
matter immersed in yS-stable nuclear matter, global charge neu- 
trality being guaranteed by a lepton background. This picture 
is obvioulsy based on the assumption that the appearance of 
the charged droplets do not significantly affect the space distri- 
bution of the leptons, i.e. that the Debye screening length Ad 
be large compared to both the typical size of the droplets and 
their separation distance (Heiselberg et al 1993 Heiselberg & 
Hjort-Jensen 20001. If this condition is not satisfied the lepton 
background is distorted in such a way as to screen electrostastic 
interactions. 

The estimates of Ad reported by Heiselberg et al (1993 1 
suggest that screening effects can be disregarded if the struc- 
tures appearing in the mixed phase of quark and nuclear matter 
have typical size and separation distance of ~ 10 fm. The re- 
sults of our calculations, that will be discussed discussed later 
in this Section, show that this appears indeed to be the case. 

The second issue deserving consideration is the stability of 
the mixed phase, i.e. the question wheter or not its energy is 
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Fig. 5. As in Fig.0] but with the EOS of quark matter obtained 
using different values of the MIT bag model parameters. Upper 
panel: B = 120 MeV/fm 3 and a s = 0.5; lower panel B = 120 
MeV/fm 3 and a s = 0.. The straight line in the lower panel is 
the double tangent resulting from Maxwell cosntruction, while 
the open circles show the extrema of the coexistence region. 

lower than the energy of the cohexisting phases of nuclear and 
quark matter. 

Formation of a spherical droplet of quark matter requires 
the energy 

& D = &c + &s ■ (39) 

where the surface contribution &s is parametrized in terms of 
the surface tension cr according to 

3 2 



o-4ttR z 



(40) 



R being the droplet radius. The electrostatic energy &c can be 
cast in the form (Ravenhall et al ll983l Pethick et al !1995> 



F 3 Q 2 (l 3 1/3 
Or = 1 U 

C 5 R\ 2 



1 

H U 

2 



(41) 



with u - (R/R c ) 3 , R c being the radius of the Wigner-Seitz cell. 
Note that the first term in the right hand side of the above equa- 
tion is the self energy of a droplet of radius R and charge Q 
obtained from Gauss law. The electric charge Q is given by 

Q = (pqm - Pnm) , (42) 

Pqm and p NM being the charge densities of quark matter and 
nuclear matter, respectively. Minimization of the energy den- 
sity € = 3£o/47i7? 3 with respect to the droplet radius yields 

&s = 2B C (43) 
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and 



R = 



4-MPqm - Pnm) 2 



3cr 



■1/3 



where 



f 3 (u) = - (2- 3k 1/3 + u) . 



(44) 



(45) 



As the density increases, the droplets begin to merge and 
give rise to structures of variable dimensionality, changing firts 
from spheres into rods {spaghetti) and eventually into slabs 
(lasagna). At larger densities the volume fraction occupied by 
quark matter exceeds 1 /2, and the role of the two phases is 
reversed. Nuclear matter, initially arranged in slabs, turns into 
rods and shperes that finally dissolve in uniform charge neutral 
quark matter. 

The energy density needed for the formation of the struc- 
tures appearing in the mixed phase has been obtained by 
Ravenhall et al ( 1983 1 in the context of a study of matter in 
the neutron star inner crust. It can be cast in the concise form 



a 2j2( NM QM\ 2 r , \ 

—cr d \n B - n~ ) f d (u) 



(46) 



where a is the fine structure constant, u is the volume fraction 
occupied by the less abundant phase (i.e. u = x f° r X < 1/2, 
u — 1 — xfoix - 1/2) and 



fd(u) = 



1 



d + 2 



d-2 



(47) 



For d = 1 , 2 and 3 Eas. (l46> and (1471 provide the correct energy 
density for the case of slabs, rods and spheres, respectively. 

For cr — both surface and Coulomb energies vanish, and 
the energy density of the mixed phase is given by Ea.(l29>. 
while for cr + 



e(cr) = e(<x = 0) + e c + e s 



(48) 



The mixed phase is energetically favorable if e(cr) is less than 
the energy density obtained from Maxwell construction, given 
by Eq.Q. 

The value of the surface tension at the interface between 
nuclear and quark matter is not known. It has been estimated 
using the MIT bag model and ignoring gluon exchange (Berger 
& Jaffe 1987 1199U . Assuming a strange quark has mass of 
~ 150 MeV Berger & Jaffe predict cr ~ 10 MeV/fm 2 . To quan- 
titatively investigate the stability of the mixed phase, we have 
calculated A e = e(cr) - e(0) for different values of cr, ranging 
from 2 MeV/fm 2 to 10 MeV/fm 2 . 

For any given value of the baryon number density «b, the 
energy density of Eas.( l46> -( I47> has been calculated using the 
nuclear and quark matter densities determined according to the 
procedure described in Section[3]and carrying out a numerical 
minimization with respect to the value of the dimensionality 
parameter d. As n B increases, the resulting values of d change 
initially from ~ 3 to ~ 2 and ~ 1 and then again to ~ 2 and 
finally to ~ 3. For example, in the case illustrated by Fig. [6] 
and corresponding to cr = 10 MeV/fm 2 , we find that sherical 
droplets of quark matter (d ~ 3) appear at ng ~ .75 firT 3 and 
turn into rods (d ~ 2) and slabs (d ~ 1) at rig ~ .95 and ~ 1.2 



frrT 3 , respectively. For larger densities quark matter becomes 
the dominant phase (i.e.;^ > 1/2): at tig ~ 1.5 and ~ 1.7 frrT 3 
the mixed phase features rods (d ~ 2) and droplets (d ~ 3) 
of nuclear matter that eventually dissolve in the quark matter 
background. 
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Fig. 6. The solid lines correspond to the difference A e = 
e(cr) - e(0) (see Eq.gH)), evaluated for cr =_10, 5 and 2 
MeV/fm 3 . The dashed line shows the difference A between the 
energy density resulting from Maxwell conatruction and e(0). 
The arrows mark the limits of the cohexistence region. Nuclear 
and quark matter are described by the APR EOS and the MIT 
bag model EOS, with a s = 0.5 and B = 200 MeV/fm 3 , respec- 
tively. 
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Fig. 7. Same as in Fig. |6] but with a s = and B = 120 
MeV/fm 3 . The solid lines show to the difference A e = e(cr) - 
e(0) (see Ea.(l48t). evaluated for cr = 5 and 2 MeV/fm 3 , respec- 
tively. 



Our results are summarized in Figs. [5] and that corre- 
spond to different choices of the MIT bag model parameters. 
The solid lines show the n B dependence of the calculated A e for 
different values of the surface tension cr. The dashed line repre- 
sents the difference A e = £m - e(0), where e« is the energy den- 
sity obtained from Maxwell construction. For any given value 
of the surface tension, the mixed pase is favorable if the corre- 
sponding solid line lies below the dashed line. 
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The results of Fig.|6j corresponding to B — 200 MeV/fm 3 
and a s = 0.5, show that the mixed phase, while being always 
the lowest energy phase for <x = 2 MeV/fm 2 , becomes to be en- 
ergetically unfavorable at some densities for <x > 5 MeV/fm 2 . 
For <x = 10 MeV/fm 2 cohexistence of charge neutral phases of 
nuclear and quark matter turn out to favorable over the whole 
density range. 

To gauge the dependence upon the MIT bag model pa- 
rameters we have repeated the calculations setting B = 120 
MeV/fm 3 and a s = 0. The results of Fig.0show that for cr in 
the range 2-5 MeV/fm 2 the mixed phase is energetically favor- 
able over a density region larger than the cohexistence region 
predicted by Maxwell construction. 

Finally, we return to the problem of the comparison be- 
tween the Debye screening length and the typical size of the 
structures appearing in the mixed phase. Our results suggest 
that the condition outlined in the work of Heiselberg et al 
( 1993 1 are indeed fulfilled. For example, in the case B = 200 
MeV/fm 3 and a s = 0.5 we find that in the region of ^ <sc 1, 
corresponding to formation of droplets of quark matter, the 
droplets radius given by Eq.J44> is ~ 2 - 3 fm. 

5. Implications for star structure 

Plugging the EOS P = P(e) into the Tolman Oppenheimer 



(49) 



maximum mass, M max ~ 1.4 M Q , barely compatible with the 
measured neutron star masses. 



Volkoff (TOV) equations (Tolman 1934 Oppenheimer & 
Volkoff tT539Tl 

dP(r) _ le(r) + P(r)l [M{r) + 4nr 2 P(r)] 

dr ~ ~° r 2 [l -2GM(r)/r] ' 

where G denotes the gravitational constant, and 



M(r) = 4n I r' z dr'e(r') 
lo 



(50) 



one can obtain the properties of the stable configurations of a 
nonrotating neutron star. Eqs.j49l> and d50i are solved by inte- 
grating outwards with the initial condition e(r = 0) = e c . For 
any given value of the central density, e c , the star radius R is 
determined by the condition P(R) = and its mass M = M(R) 
is given by Eq.dBol. 

The occurrence of the transition to quark matter makes the 
EOS softer, thus leading to a lower value of the maximum 
mass. In Fig. [8] we compare the mass-central energy density 
relations obtained using the APR EOS only to that obtained al- 
lowing for a transition to quark matter with a s = 0.5 and B = 
120 and 200 MeV/fm 3 . The transition is described according to 
Gibbs conditions, neglecting surface and Coulomb effects. We 
find M max = 2.20 M Q for the star made of nuclear matter oly 
and M max = 1.89 and 2.03 M for the hybrid stars correspond- 
ing to B = 120 and 200 MeV/fm 3 , respectively. 

In Fig. [S] we also compare the M{e c ) curves obtained set- 
ting B — 120 MeV/fm 3 and a s — and adopting either 
Gibbs or Maxwell picture. Whether the phase transition pro- 
ceeds through the appearance of a mixed phase or charge neu- 
tral cohexisting phases does not appear to significanltly affect 
the mass-central energy density relation. On the other hand, 
neglecting perturbative gluon exchange results in a rather low 




Fig. 8. Relation between neutron star mass and central energy 
density for different EOS. Dotdash line: pure nuclear matter 
(APR EOS); dotted lines: nuclear matter (APR EOS) and quark 
matter (MIT bag model with B = 120 and 200 MeV/fm 3 and 
a s = 0.5). The phase transition is described according to Gibbs 
rules; solid line: same as the dotted line, but with B = 120 
MeV/fm 3 and a s = 0.5; dashed line: same as the solid line, but 
with the phase transition described using Maxwell construc- 
tion. 



The neutron star mass-radius relations associated with the 
M(e c ) curves of Fig. |8j displayed in Fig. |9j show that in this 
case using Maxwell construction instead of Gibbs rules pro- 
duces a visible effect. It has to be noted that all EOS predict 
the existence of stable star configurations with masses in the 
range allowed by observation (Thorsett & Chakrabarty 1999 
Qaintrell et al 2003}, as well as a M(R) relation compatible 
with that resulting from the gravitational red shift measurement 
ofCottam et al 
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Fig. 9. Mass radius relation for different EOS. The meaning of 
the curves is the same as in Fig. [8] The horizontal lines corre- 
spond to the observational limits on neutron star mass, whereas 
the third straight line is the mass-radius relation resulting from 
the gravitational redshift measurement of Cottam et al (2002i. 
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Fig[IT"]shows that different descriptions of the phase transi- 
tion lead to remarkably different star density profile. While in 
presence of the mixed phase the density is a smooth function of 
the distance from the star center, Maxwell construction leads to 
the appearance of a disontinuity. For comparison, we also show 
the profile of a star of the same mass, ~ 1 .4 M , made of pure 
nuclear matter described by the APR EOS. 

The discontinuous behavior can be easily understood not- 
ing that TOV equations i49i and ( 15 01 require that the pressure 
P(r) be a monotonically decreasing function. It follows that if 
the pressure is the same for two different values of density, as 
in the phase transition a la Maxwell, they must necessarily cor- 
respond to the same value of r. 




2 4 6 8 10 12 

r [Km] 



Fig. 10. Density profiles of a neutron star of mass ~ 1 .4 M re- 
sulting from different models. Dashed line: APR EOS and MIT 
bag model EOS with a s = and B = 120 MeV/fm 3 . Phase tran- 
sition described allowing for the presence of the mixed phase. 
Solid line: same as the solid line, but with the phase transition 
described according to Maxwell construction. Dotdash line: 
pure nuclear matter with APR EOS. 

In the cohexisting phases scenario, the transition only takes 
place in stars whose central density exceeds the density of the 
quark matter phase. These star configurations turn out to be 
marginally stable, their mass being close to the maximum mass. 
For example, setting B = 200 MeV/fm 3 and a s = 0.5 we find 
that the transition only occurs in stars having mass ~ 2.OM . 
The radius of the quark matter core is pretty small (~ 1 Km), 
while the density jump is large, going from 1.8 x 10 15 g/cm 3 to 
3.9 x 10 15 g/cm 3 . These results are to be compared to those ob- 
tained in the mixed phase scenario when Coulomb and surface 
effects are neglected. In this case there is no jump and the den- 
sity varies smoothly. At the center of a star of mass ~ 2.OM , 
corresponding to energy-density ~ 2.4 x 10 15 g/cm 3 the volume 
fraction occupied by quark matter reaches ;f ~ 30 %. 

6. Conclusions 

We have carried out a study of the transition from nuclear mat- 
ter to quark matter in the inner core of neutron stars, aimed at 
assessing whether the appearance of a mixed phase is energeti- 
cally favorable and how the emerging picture depends upon the 
parameters entering the MIT bag model EOS. 



In order to minimize the uncertainty associated with the de- 
scription of the nuclear matter phase, we have adopted a EOS 
obtained from an ab initio calculation, based on a dynamical 
model stongly contrained by experimental data and not involv- 
ing any adjustable parameters (Akmal et al 1998 1. 

Our results show that the effect of perturbative gluon ex- 
change on the MIT bag model EOS is large and cannot 
be accounted for adjusting the value of the bag constant B. 
Neglecting interactions between quarks leads to a considerable 
softening of the EOS, resulting in a drastically lower transition 
density. 

The fact that the setting a s = lead to a much softer EOS 
is reflected by the rather small value of the maximum neutron 
star mass, barely exceeding the canonical value of 1.4 M . 

The stability of the mixed phase turns out to be strongly af- 
fected by surface and Coulomb effects. Using the softer quark 
matter EOS (a s - 0,B - 120 MeV/fm 3 ) we find that even a 
very small value of the surface tension, <x = 2 MeV/fm 3 , pro- 
duces a narrowing of the density region spanned by the mixed 
phase. With the harder EOS (a s = 0.5, B = 200 MeV/fm 3 ) the 
cohexistence of neutral phases of nuclear and quark matter is 
energetically favored at all densities for cr — 10 MeV/fm 3 . 

Comparison with the results of Heiselberg et al suggests 
that surface and Coulomb effect become larger when a realistic 
EOS is employed to describe the nuclear matter phase. 

While the mass-central density relation appears to be 
largely unaffected by the occurrence of the mixed phase, its 
effect can be clearly seen in the M(R) curve. However, as mea- 
surements of neutron star radii are plagued by large uncertain- 
ties, this feature is not likely to be exploitable to extract clearcut 
information from observations. 

The most striking difference between the Maxwell and 
Gibbs picture of the phase transition shows up in the neutron 
star density profile, which in the case of transition at constant 
pressure exhibits a sharp discontinuity. 

The presence of a density jump is known to affect neutron 
star dynamics, leading to appearance of a class of nonradial 
oscillation modes, called g-modes, associated with emission of 
gravitational radiation. 

Early investigations of the g-modes focused on the discon- 
tinuities produced by the changes of chemical composition in 
the low density region of the neutron star crust, corresponding 
to a fractional distance from the surface < 10 % (Finn 1198*71 
McDermott 1990 1. These studies have been recently extended 
to the case of g-modes produced by a discontinuity located at 
much larger density and involving a much larger density jump, 
such as those associated with the transition to quark matter 
(Miniutti et al 2003 1. Based on the results of calculations car- 
ried out using siple polytropic EOS, Miniutti et al (2003 1 argue 
that a simultanuoeus measurements of the frequencies of the 
fundamental /-mode and the g-mode would provide informa- 
tion on both size and location of the discontinuity. 

Although it is unlikely that the first generation of laser in- 
terferometric antennas will detect gravitational waves emitted 
by an oscillating neutron star, the new detectors currenltly un- 
der investigation (see, e.g., the EURO proposal (2000$) are ex- 
pected to be much more sensitive at the relevant frequencies 
above 1-2 kHz. Hopefully, information on the neutron star mat- 
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ter EOS and the nature of the transition to quark matter may 
come gravitational wave astronomy. 
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